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ABSTRACT 

We carry out classification of 4330 X-ray sources in the 2XMMLDR3 catalog. They are selected 
under the requirement of being a point source with multiple XMM-Newton observations and at least 
one detection with the signal-to-noise ratio larger than 20. For about one third of them we are able 
to obtain reliable source types from the literature. They mostly correspond to various types of stars 
(611), active galactic nuclei (AGN, 753) and compact object systems (138) containing white dwarfs, 
neutron stars, and stellar-mass black holes. We find that about 99% of stars can be separated from 
other source types based on their low X-ray-to-IR flux ratios and frequent X-ray flares. AGN have 
remarkably similar X-ray spectra, with the power-law photon index centered around 1.91±0.31, and 
their 0.2-4.5 keV flux long-term variation factors have a median of 1.48 and 98.5% less than 10. In 
contrast, 70% of compact object systems can be very soft or hard, highly variable in X-rays, and/or 
have very large X-ray-to-IR flux ratios, separating them from AGN. Using these results, we derive 
a source type classification scheme to classify the other sources and find 644 candidate stars, 1376 
candidate AGN and 202 candidate compact object systems, whose false identification probabilities 
are estimated to be about 1%, 3% and 18%, respectively. There are still 320 associated with nearby 
galaxies and 151 in the Galactic plane, which we expect to be mostly compact object systems or 
background AGN. We also have 100 candidate ultra-luminous X-ray sources. They are found to be 
much less variable than other accreting compact objects. 
Subject headings: catalogs — X-rays: general — infrared: general 



1. INTRODUCTION 

The XMM-Newton observatory has been carrying out 
pointed observations of the X-ray sky for more than a 
deca de. The XMM-New ton Serendipitous Source Cata- 
log (|Watson et alj [20091 ) provides various properties of 
serendipitously detected X-ray sources, in addition to 
targets, from these pointed observations. It has been 
updated regularly, and the latest version as of the writ- 
ing of this paper is the 2XMMLDR3 catalog, containing 
353191 X-ray source detections for 262902 unique X-ray 
sources. It is the largest X-ray source catalog ever pro- 
duced. 

Such catalogs are rich resources for exploring a va- 
riety of X-ray sour ce populations and identifying rare 
source types (e.g . . iLiii et al.1 120111: jPineau et al.l 120111 : 
iPires et all l2009UFarrell et al.l 120091) . Compared with 
other large X-ray catalogs such as the ROSAT cata- 
logs from its all-sky survey (RASS) and p ointed obser- 
vatio n s and the Chandra Source Catalog (iWhite et al.1 
fl99"i IVoges et all [l999l [200(1 : lEvans et al.ll2010D . the 
2XMMLDR3 catalog has several advantages because of 
the large effective area, high spatial resolution, and/or 
broad energy band coverage of the XMM-Newton obser- 
vatory. 

In this paper, we select a sample of sources of the 
best quality from the 2XMMLDR3 catalog and study 
their properties in detail. We first establish a subsam- 
ple of these sources whose types are known and study 
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their X-ray spectral shapes, X-ray variability, and X-ray- 
to-optical and X-ray-to-IR flux ratios. Based on these 
properties, we define some quantitative criteria to clas- 
sify other sources. Sources with interesting properties, 
discovered as a part of this work are the subject of a 
companion paper. We devote significant effort to visu- 
ally screen the results so as to reduce various kinds of 
systematic errors. 

In Section [5J we describe the source selection proce- 
dure, the calculation of the flux and the hardness ratio, 
the search for stellar X-ray flares, the measurement of the 
long-term variability, the source classification from the 
literature, the search for optical and IR counterparts, the 
determination of the (extra) galactic nature, and simple 
spectral fits. The properites of the sources with known 
types are shown in Section 13.11 The scheme that we 
propose to classify sources is given in Section 13.21 The 
results of applying this scheme to the sources with un- 
known types are presented in Section ^. 31 The discussion 
and conclusions of our study are given in Section |4j 

2. DATA ANALYSIS 

2.1. The 2XMMi-DR3 Catalog and Source Selection 

The 2XMMi-DR3 catalog is based on 4953 pointed ob- 
servations. These observations have a variety of modes 
of data acquisition, exposures (>1000 s), and optical 
blocking filters (thick, medium, and thin). The cata- 
log contains detections in five basic individual energy 
bands (0.2-0.5, 0.5-1.0, 1.0-2.0, 2.0-4.5, and 4.5-12.0 
keV, numbered 1-5, respectively), for each European 
Photon Imaging Camera (EPIC) i .e., pn, MPS 1 /Ml 
and MOS2/M2 (IJansen et^[2TOlt Intruder et al.M200Lb 
iTurner et all I2001D . Results are also given for various 
combinations of energy bands and cameras. Band 8 
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refers to the total band 0.2-12.0 keV, and EPIC/EP 
refers to all cameras. We also follow these notations 
throughout this paper. 

We focused only on point sources with multiple XMM- 
Newton observations and at least one detection with the 
signal-to-noise ratio S/N >20 in this paper. We consid- 
ered only point sources because extended sources (typ- 
ically Galactic supernova remnants (SNRs) and galaxy 
clusters) tend to be subject to complicate d problems of 
sourc e detection and flux measurement (jWatson et al.l 
2009). Wc dchnc a source to be point-like when its 
extent (EP_EXTENT in the catalog) is zero (at least in 
one detection). We note that in the pipeline, if the ex- 
tent radius in a detection is <6", it is set to zero, i.e., 
a point source is assumed. We estimated the S/N us- 
ing EP_8_CTS/EP_8_CTS_ERR, where EP_8_CTS, given in 
the catalog, is the total counts summed over all cam- 
eras, and EP_8_CTS_ERR is its error. Our requirement of 
having at least one detection with S/N >20 is to en- 
sure at least one good spectrum for detailed analysis 
when needed. The requirement of having multiple ob- 
servations allows us to calculate the long-term variabil- 
ity. It is different from the requirement of having mul- 
tiple detections. The catalog includes only detections 
with a total likelihood (from all bands and all cameras) 
> 6 (corresponding to about 3a). We used the service 
http://www. ledas.ac.uk/flix/flix. html to find the 
observations of a source and estimate the fluxes for the 
observations that have no entries in the catalog. 

There are 5802 sources satisfying the above conditions. 
However, as the catalog was produced using an auto- 
mated procedure, some of these sources suffer from var- 
ious kinds of problems, and we excluded them based on 
visual inspection. About 16% are spurious, mostly due 
to single reflections caused by bright sources outside the 
field of view, out-of-time ey ents, and extremely bright 
sources (| Watson et al.l l2009h . We also excluded about 
8% of the sources that are detected within bright ex- 
tended sources, especially SNRs, though they pass the 
above criteria. About 2% of the sources are finally ex- 
cluded due to problems such as optical loading, serious 
event pileup, and being too close to other bright sources. 

In the catalog, each source has a unique source number 
(SRCID in the catalog). However, from visual inspection, 
we still found 42 sources assigned with multiple SRCID 
numbers (Table [l}. Each source will be referred to using 
the first number in Table [TJ but its detections consist of 
all the detections corresponding to the different SRCID 
numbers (Table [3]). In the end, we have 4330 unique 
sources in our sample (Table 2]). 

2.2. Calculation of the Flux and Hardness Ratio 

The 2XMMi-DR3 catalog provides the count rates and 
fluxes of each detection for various combinations (indi- 
vidual and total bands/cameras). The count rates were 
corrected for various instrumental effects including the 
mirror vignetting, detector efficiency and point spread 
function losses. The fluxes in individual energy bands 
were calculated by dividing the count rates by the en- 
ergy conversion factors (ECFs). The ECFs depend on 
the camera, the filter, and the source spectrum, which 
the catalog assumed to be an absorbed power-law (PL) 
with a photon index of Tpl = 1-7 and an absorption 
column density of TVh = 3 x 10 20 cm~ 2 (jWatson et all 



TABLE 1 

Sources with multiple SRCID numbers 



106501 241552 


111034 111032 


113853 113854 


127424 127426 


134832 134833 210662 


13940 195048 


143291 143295 


147800 147801 


150531 150532 


163538 163539 


183730 183731 


190860 190861 


198805 42983 42980 


201208 201209 


202064 62202 


20302 196589 


205954 86092 


206123 206125 


20843 196637 


210401 245122 


213162 149986 


221841 4783 


227999 198991 


235153 73445 


235311 73852 


252232 156258 


37588 37584 


49466 49454 


50074 50083 


50260 50265 


50339 50335 


53515 53514 


54078 54077 


61618 61621 


6320 6324 


63368 63365 


66168 66171 


69639 69640 


73642 73647 


8091 8087 


91437 206280 206281 


93326 93329 





Note. — Each cell lists the SRCID numbers that should corre- 



spond to the same source. 

2009). The ECFs used in the catalog were based on the 
calibration before 2008, except for a few observations. 
We applied correction factors provided on the catalog 
website to obtain fluxes corresponding to the newer cal- 
ibration. In this paper, we mostly used the EPIC fluxes 
in 0.2-12.0 keV (band 8) and 0.2-4.5 keV (called band 
14 hereafter). They were obtained by first summing the 
fluxes in the relevant individual bands and then averag- 
ing over the available cameras weighted by the errors. We 
estimated their systematic errors to be 0.107 and 0.074, 
respectively (Appendix [XJ . 

We also calculated four X-ray hardness ratios HR1 
HR4 defined as 

URi = {R l+1 - Ri)/(Ri+i + Ri), (1) 

where Ri and Ri+i are the MOSl-medium-filtcr equiva- 
lent count rates in the energy bands i and i + 1, respec- 
tively. These count rates were obtained by multiplying 
the EPIC fluxes by the ECFs for the MOS1 camera with 
the medium filter. We note that the 2XMMi-DR3 catalog 
also provides the camera-specific X-ray hardness ratios, 
but the count rates used, and thus the hardness ratios, 
depend on the camera and the filter. This poses a prob- 
lem to use the hardness ratios to investigate the source 
spectral shape uniformly, considering that the observa- 
tions used to create the catalog have various filters and 
available cameras. Our hardness ratios are less subject 
to this problem, as the fluxes have been calculated to 
mitigate their dependence on the camera and filter by 
using the camera- and filter-dependent ECFs. We only 
observed small systematic differences between the pn and 
MOS hardness ratios using our definition (Appendix [A"|. 

2.3. Variability Calculation 

The variability properties can provide important clues 
about the source type. There are many kinds of variabil- 
ity. For the short-term variability that can be measured 
within a single observation, we paid attention to stellar 
X-ray flar es, which are good indicators of coronally ac- 
tive stars (|Gudelll2004f ). They show a variety of profiles, 
though most of them have a fast rise and a slow decay 
and can last from minutes to hours. The light curves 
created for bright detections by the pipeline are suitable 
for stellar X-ray flare search. They are extracted from a 
circular region of radius 28" and have bin sizes an integer 
multiple of 10 s ( the minimum is 10 s ), depending on the 
source intensity (jWatson et al.l [2009). We describe our 
search for stellar X-ray flares in Appendix [Bj 
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TABLE 2 
Source type statistical breakdown 



Type 


# 


Selection 


Main Reference 


Identified: 


Star 


202 






OrVr 


252 






PrSt 


27 




SIMBAD, GCVS 


VrSt 


101 






FISt 


29 






Syln 


29 


Sp="Sln" 




Syl 


242 


Sp="Sl" 




Sy2 


63 


Sp="S2", "Sih" or "Sli" 




LIN 


12 


Sp="S3", "S3b" or "S3h" 


VV10 


Bla 


27 


Cl= "B" 




QSO a 


250 


01= "Q" 




AGN a 


130 


Cl= "A" except Sp="H2" 




INS 


7 


Typc="XDINS" 




MGR. 


11 


Type="AXP" 


ATNFPC 


rPsr 


21 


Others 




aPsr 


34 






Bstr 


18 




Literature 


BHB 


4 






CV 


43 




CCB 


Candidate: 


Star 


644 




This work 


AGN 


1292 




This work 


G 


84 




This work 


SNR 


16 




Literature 


Mixed 


19 




Literature 


ULX 


100 




Literature 


CO 


202 






XGS 


320 




This Work 


GPS 


151 






Note. 


— See 


Scctionsl2.4l and 13.31 for the detailed description 



of each source type. 

a Excluding identified Scyfcrt galaxies, LINERs 



bla 



We also measured the long-term flux variability over 
different observations. Both bands f4 and 8 were used, 
but we will explore results with band f4 in more detail 
because the flux in band 5 (4.5-12.0 keV) tends to have 
large uncertainties. We defined the long-term flux vari- 



ability as V V i 
difference as 

S v. 



-Fmax/-Fmin and the significance of the 



F„ 



2 

max 



{rF r 



) 2 

max ) 



(rF a 



\2\l/2' 



(2) 



EPIC fluxes of a unique source, with the corresponding 
(statistical) errors cr ma x and a m - ul , respectively, and r is 
the systematic error (r = 0.107 and 0.074 for bands 8 and 
14, respectively; Section |2~21 . We only used detections 
with the flux above four times the error (cr) when calcu- 
lating -Fmax, while we used 2a as the flux for detections 
with the flux less than 2a when calculating i^ m i n - 

2.4. Source Type Identification from the Literature 

We established a subsample of our source list with well 
confirmed source types from the literature (Tables [2] and 
as described below. We refer to them as identified 
sources hereafter. There are three main categories of 
X-ray sources: stars, active galactic nuclei (AGN) and 
compact object systems containing white dwarfs (WDs), 
neutron stars (NSs) and stellar-mass black holes (BHs). 

For stars, we used the General Catalog of Variable 
Stars (GCVS, Version 2011 January. ISainus et al.ll2009l) 
and sources from SIMBAD (three from the literature) 
with optical spectral types available. We classified stars 



into the following categories: Orion variables (OrVr); 
pre-main-sequence stars (PrSt); variable stars (VrSt, ex- 
cluding Orion variables and flaring stars, mostly eclips- 
ing binaries, rotationally variable and pulsating variable 
stars); flaring stars (FISt); and stars for the rest, typi- 
cally main-sequence stars. 

For A GN, we used the catalog of q uasars and active 
nuclei by iVeron-Cettv k Veronl pM VV10 hereafter), 
with an additional three from the literature (based on 
optical spectra) . We classified them into seven main cat- 
egories based on the classifications given in this cata- 
log (Table [2]): narrow-line Seyfert 1 (Syln); Seyfert 1 
(Syl); Seyfert 2 (Sy2); low-ionization nuclear emission- 
line regions (LINERs or LINs); blazars (Bla); quasi- 
stellar objects (QSOs, i.e., quasars); and AGN for the 
rest. There are 11 Seyfert galaxies with broad polarized 
Balmer lines (Sp="Slh" in VV10) or broad Paschen lines 
in the infrared (Sp="Sli"). They are generally classified 
as Seyfert 2 galaxies in SIMBAD, and we adopt the same 
scheme. 

The Australia Telescope National F acility Pulsar Cat- 
alog (ATNFPC, as of 2011 March, Manchester et al.l 
20050 was used to identify thermally cooling isolated 
NSs (INSs), rotation-powered pulsars (rPsr), and mag- 
netars (MGRs). The accretion-powered X-ray pulsars 
(aPsr) and bursters (Bstr, believed to be weakly magne- 
tized accreting NSs (mostly low- mass X-ray binaries)) 
were identified from the literature, with the require- 
ments of detections of pulsation s and type-I X-ray bursts, 
respec tively. We referred to [McClinto ck k Remillardl 
(2006) to identify four BH X-ray binaries (BHBs), which 
either have secure BH mass measurements or show 
X-ray properties very similar to those of BHBs with 
secure mass measurements (g rade A in Table 4.3 in 
iMcClintock fe Remillardl 120061) . The Catalog of C ata- 



clysmic Binaries (CCB) from lRitter k Kolbl poM 7th 
edition) was used to identify the cataclysmic variables 
(CVs). All these types of sources will be generally re- 
ferred to as compact objects (COs) in this paper. 

We also found 16 sources that are probably (Galactic or 
extragalactic) SNRs (or their knots) from the literature. 
The 19 "mixed" sources in Table [2] (see also Table HJ) 
include two cooling WDs, three supernova, four gamma- 
ray bursts, four symbiotic stars, two micro-quasars, etc. 
The 100 ultra-luminous X-ray sources (ULXs, with lumi- 
nosity above 10 39 erg s _1 ) are either from the literature 
or from this work (Section I3.3|) . We treated all these 
SNRs, "mixed" objects, and ULXs as candidate (instead 
of identified) sources, considering that their X-ray mech- 
anisms are largely still unclear. 

2.5. Multi-wavelength Cross-correlation 

We cross-correlated our source X -ray p ositions with the 
USNO-B1.0 Catalog (I Monet et al. l l2T)03h and the 2 MASS 
Point Source Catalog (2MASS PSC. ICutri eTaII[2003l) 
to search for their optical and IR counterparts, respec- 
tively. Our sources have small positional errors, with 
86% less than 0"5 and 99% less than l'/O (statistical plus 
systematic). The USNO-B1.0 Catalog and the 2MASS 
PSC have comparable positional errors. We selected the 
counterpart to be the closest one within 4". For a few 

4 http:/ /www. atnf.csiro.au/research/pulsar/psrcat 
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sources, we did not assign the counterparts in this way. 
One common case is bright extended galaxies, for which 
the USNO-B1.0 Catalog and the 2MASS PSC often have 
many entries, and we chose the brightest nearby one (in 
terms of the i?2/i^ s -band magnitude) as the counterpart. 
The other common case is stars with large proper mo- 
tions (indicated in the USNO-B1.0 Catalog). In the end 
we found optical and IR counterparts for 3014 (70%) and 
2058 (48%) of our sources, respectively (Table H}. For 
very few sources, we denoted the match quality as bad 
due to their coincidence with the persistence artifact po- 
sitions in the 2MASS0, bright globular clusters, etc. and 
we will not explore their optical/IR properties. 

We used the i?2-band magnitude to define the op- 
tical flux as logpo) = -B2/2.5 - 5.37, following 
iMaccacaro et alJ (fr988). In some cases when the R2- 
band magnitude is not available, we used the i?l-band 
magnitude instead. For sources without optical counter- 
parts, we assumed R2 — 21.0, which is about the min- 
imum among the counterparts found and was taken as 
an upper limit. The IR flux was calculated with the 
ifs-band magnitude from the 2MASS PSC: log(F m ) = 
—K s /2.5 — 6.95, assumi ng the "in-band" ze ro-magnitude 
flux of the K s band from lCohen e t al. (2003|) . For sources 
without IR counterparts, we calculated the upper limit 
using K s = 15.3, the 3-er limiting sensitivity of this 
band in the 2MASS. Throughout the paper, all fluxes 
are in units of erg s _1 cm -2 and correspond to appar- 
ent/absorbed values. 

We also similarly cross-correlated our sources with the 
extended IR sources fr om the 2MASS Exte nded Source 
Catalog (2MASS XSC. ISkrutskie et alJ2006h and the op- 
tical sources from the Sloan Digital Sky Survey (SDSS, 
lAbazajian et al.ll2009f ). We used them mostly to check 
whether our sources are coincident with the center of (ex- 
tended) galaxies. To reduce the contamination, we only 
assumed extended IR sources from the 2MASS XSC to 
be galaxies when they are coincident with the center of 
an optical extended source, based on visual inspection of 
the Digital Sky Survey (DSS) images. For the SDSS ex- 
tended sources, we assumed them to be galaxies if their 
r-band magnitudes are >20 (the separation between the 
point and extended sources by t he SDSS pipeline is n ot 
reliable for very faint sources (see lScranton et aT] |2002)). 

2.6. Association with Nearby Galaxies 

We used the Third Reference Catalog (RC3) of galax- 
ies (|de Vaucouleurs et alj|199il ) to investigate the prob- 
ability of our sources being non-nuclear extragalactic 
sources. This catalog is reasonably complete for galax- 
ies having apparent diameters larger than V at the D25 
isophotal level and total B-band magnitudes brighter 
than about 15.5 and redshifts less than 0.05. The -D25 
isophote refers to the 25 mag arcsec -2 i?-band isophote 
and is approximated as an ellipse. This catalog also in- 
cludes some objects of special interest, resulting in 23,022 
galaxies in total. It gives the apparent major and minor 
diameters and the position angle of the -D25 isophote for 
each galaxy. The D25 isophote is roughly the domain 
of a galaxy as seen in the DSS images, and by compar- 
ing the source positions with this isophote we can check 
whether the sources are possibly associated with that 



galaxy. Following iLiu fc Bregmanl (|2005[) , we calculated 
the angular separation a between the galaxy center and 
the source, which was then scaled by the elliptical radius 
i?25 of the D25 isophotal ellipse in the direction from 
the galaxy center to the source. We assumed that only 
sources with a/i?25<2 and a>4" can possibly be non- 
nuclear extragalactic sources. The positions of the galaxy 
centers given in the RC3 have only modest accuracy. We 
found that the positions from NASA/IPAC Extragalac- 
tic Database (NED)0 are closer to the real galaxy cen- 
ters (at least for galaxies related to our sources) based 
on visual inspection of DSS images. Thus we chose to 
use the galaxy center positions from NED when possible. 
We also estimated the maximum 0.2-12.0 keV absorbed 
luminosity for each candidate non-nuclear extragalactic 
so urce. We adopted the d ista nces t o mos t galaxies used 
in ILiu fc Bregmanl ((20051 ) and Eh] (|20ll . For 17 other 
galaxies, we obtained their distances using the redsifts 
from NED and assuming a flat universe with the Hubble 
constant Ho=75 km s _1 Mpc -1 and the matter density 
f-M=0.27 if they have recessional velocities larger than 
1000 km s _1 or from the literature otherwise. 

2.7. Simple Spectral Fits 

We created simple energy spectra, one for each avail- 
able camera for each detection, using the count rates in 
five basic energy bands in the 2XMMi-DR3 catalog and 
fitted them with an absorbed PL (Table [3]). The re- 
sponse files corresponding to the "thin" , "medium" , and 
"thick" optical blocking filters were constructed using the 
observations 0112560101, 0204870101, and 0311190101, 
respectively, from source regions centering at the point- 
ing direction and covering the whole field of view. SAS 
11.0.0 and the calibration files of 2011 June were used. 
The event selection criteria for each energy ban d / camera 
used for source detection in lWatson et al I (|2009ft were fol- 
lowed. The interstellar medium absorption was modeled 
by the WABS model in XSPEC. A systematic error of 
10% on the model was assumed. Relative normalizations 
among different cameras were allowed to be free. We note 
that we used the PL model just to roughly characterize 
the spectral shape. The real detailed X-ray spectra of 
AGN are still often described by this model when there 
are no strong soft excesses, but for stars, the APEC and 
MEKAL models in XSPEC are generally needed to de- 
scribe their detailed X-ray spectra well. 

3. RESULTS 

3.1. Properties of Identified Sources 

3.1.1. Spatial distribution 

We show the spatial distribution of the identified 
sources in Galactic coordinates in Figure [1] (left panel) 
and their distributions with respect to the Galactic lat- 
itude b (right panel). We differentiate stars, AGN, and 
compact objects. We note that in the sky map in the left 
panel some fields are too crowded to see the real num- 
ber of sources, while the distributions on the right panel 
complement this information. The bin sizes in the distri- 
bution plots vary to correspond to an equal spatial area. 
The error bars shown were obtained by assuming Poisson 
statistics in each bin (this is also assumed for all other 



5 http:/ /www. ipac.caltech.edu/2mass/releases/allsky/doc/sec4_7. html 6 http:/ /ned. ipac.caltech.edu 
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Fig. 1. — Hammer-Aitoff equal area projection in Galactic coordinates of identified sources (left) and their distributions with respect to 
the Galactic latitude (right). 
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Fig. 2. — The X-ray color-color diagrams for identified sources. Overplotted are PL spectra (dotted lines) with Tpl = 0.5 (top), 1, 2, 
3, and 4, and thermal spectra (APEC model with 0.5 solar abundance, dashed lines) with temperatures of 0.5 (bottom), 1, 2, and 5 keV. 
They are obtained with Nn varying from (lower-left) to 10 23 cm -2 . The detections for each source are connected by solid lines in an 
increasing order of HR3. 
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distribution plots throughout this paper). We see that 
there are many stars near b — —19.4°, which is due to 
the Orion Nebula (about 42% of the identified stars are 
from this nebula). Apart from this field, the identified 
stars concentrate at low Galactic latitudes. Stars with 
high Galactic latitudes are also seen (about 31% with 
\b\ > 20°, excluding those from the Orion Nebula). Stars 
are typically in the Galactic plane. The large spread of 
stars along the Galactic latitude in our sample indicates 
that they are nearby. 

The identified AGN tend to be at high Galactic lati- 
tudes, only about 2% with |&| < 20°. Such a bias can 
be explained by the large absorption along the Galac- 
tic plane and the preference to target the region outside 
the Galactic plane in galaxy surveys, such as the SDSS. 
There are more identified AGN from the northern Galac- 
tic hemisphere (about 68% with b > 20° while about 30% 
with b < —20°), which can be explained by the fact that 
many AGN in our sample are discovered by the SDSS, 
whose sky access is primarily in the north. 

The identified compact objects clearly concentrate 
near the low Galactic latitude region, about 47% with 
\b\ < 10°. The excess of the identified compact objects 
in the region of —50° < |&| < —30° is due to the presence 
of LMC and SMC. 

3.1.2. X-ray color-color diagrams and spectral fits 

Figure [5] shows the X-ray color-color diagrams for 
identified sources. Stars, AGN, and compact objects 
are plotted in the left, middle, and right columns, re- 
spectively. Only data points with error bars less than 
0.1 on both colors are plotted. Spectral tracks of ab- 
sorbed PL (dotted lines) and APEC spectra (dashed 
lines) are overlaid for comparison. APEC is a model 
of the emission spectrum from collisionally-ionized dif- 
fuse gas in XSPEC and is often used to model the X- 
ray spectra of stars. In the HR1-HR2 diagram, we 
see that most stars occupy a region with character- 
istic APEC temperatures about 0.5-1 keV. However, 
in the HR2-HR3 and HR3-HR4 diagrams, there are 
many stars in a region with APEC temperatures >1 
keV, which means that these stars probably need mul- 
tiple APEC components to describe them. Investiga- 
tions of these sources show that a majority of them 
correspond to detections with flares. Stars rarely show 
very hard X-ray spectra. Be stars such as HD 110432 
(#110586) and SS 397 (#1670 02) can be very hard (e.g., 
lLopes de Oliveira etafl 12007ft . with HR4 larger than 
—0.4 in Figure [5J The outliers in the upper-left region 
in the HR3-HR2 diagram consist of famous stellar sys- 
tems such as r) Carinae (#87820), j 2 Velorum (#65891), 
HL TAU (#38943 ), EX Lupi (#248785) and FU Ono- 
nis (#53697) (e g..lFMn"ier et al.ll201 ll iSchild et al.ll200l 
Giardino et al.l 120061 : IGrosso et al.1 120101 : iSkinner et al.1 
2006). 

Figure [5] shows that AGN (middle column) are mostly 
in the region with Tpl =1-3. The main outliers are 
Seyfert 2, which have HR4 larger than 0. They may 
have strong absorption and typically have two main com- 
ponents with the hard one probably due to the reflection 
and/or heavy absorption of Seyfert 1 nuclei dTurner et al.l 
[19971 lAwaki et al.ll2006l: iNoguchi et al.ll2009ft . Compar- 
ing the HR1-HR2 and HR3-HR4 diagrams, we see that 
most narrow-line Seyfert 1 galaxies are in a region with 



TpL >2 in the HR1-HR2 diagram, but in a region with 
TpL <2 in the HR3-HR4 diagram. This can be explained 
by the frequent presence of soft excesses in these sources. 

Compact objects show a large spread in Figure [2] (right 
column). This is to a large degree due to many subclasses 
included. However, many types of compact objects are 
known to show large spectral variations. For example, 
some CVs (brown stars) can have spectra changing from 
very hard to super-soft. Some weakly magnetized accret- 
ing NSs (bursters) also have a large range of HR3 (green 
triangles). Accretion-powered X-ra y pulsars, whose X - 
ray spectra are known to be hard ([White et al.l Il983h . 
occupy the region with Tpl around 0.5, where we see 
few stars and AGN. The thermally cooling isolated NSs 
(black filled circles) stay in the lower left corner of the 
(HR1-HR2 and HR2-HR3) diagrams, as expected due to 
their soft nature ([Haberll [2007'). The magnetars (pur- 
ple diamonds) look very hard in the HR1-HR2 diagram, 
which can be explained by large absorption, but many 
of them occupy a soft region in the HR3-HR4 diagram, 
with TpL >3. Many magnetars are known to show strong 
spectral turnover at 15 keV, with very soft and very 
hard spectra below and above this en ergy, respectively 
([Enoto et al.ll2010t iKuiper et al.ll2006f ). and our results 
are consistent with this. 

The spectral differences among various classes of 
sources can also be seen from the simple fits with an 
absorbed PL to spectra created using the band count 
rates in the 2XMMi-DR3 catalog (Table EJ. Focusing on 
detections with S/N > 14, we obtained 1374 (79%), 2271 
(92%) and 364 (65%) detections with reduced \ 2 < 4 for 
stars, AGN, and compact objects, respectively. Their 
TpL distributions are shown in Figure [3[ and we see that 
they are remarkably different for different classes. AGN 
have TpL concentrated near the median of 1.91, with a 
standard deviation of only 0.31, while stars and com- 
pact objects have Tpl much more scattered. The Tpl 
values for stars are generally high, with only about 13% 
less than 2.5. We found that a majority of detections 
with low TpL are due to flares, while detections with 
TpL > 6.0 (corresponding to low APEC temperatures in 
the color-color diagrams) are mostly (about 70%) from 
main sequence stars. For compact objects, detections 
with TpL > 6.0 are mostly from thermally cooling iso- 
lated NSs and novae (a subclass of CV). There are 31% 
of detections with Tpl < 1.0, mostly from accretion- 
powered X-ray pulsars. In comparison, AGN and stars 
have only 1.7% and 0.5% of detections with Tpl < 1.0, 
respectively. 

3.1.3. Long-term X-ray variability 

Figure [4] plots the X-ray flux variation factor V val ver- 
sus the maximum flux in 0.2-4.5 keV and the distribu- 
tions of V var for the three main classes (i.e., stars, AGN, 
and compact objects). We obtained median variation 
factors m — 1.79, 1.48, and 1.65 for stars, AGN, and 
compact objects, respectively. There are about 95.4% of 
the stars, 98.4% of the AGN, and 71.7% of the compact 
objects with V vaY less than 10. The distributions of V var 
for stars and AGN are relatively smooth. They roughly 
resemble half-normal distributions for the logarithm of 
Vvar (dotted lines, with the same median as the data). 
When the flux logarithm for each source follows the same 
normal distribution and there are only two detections for 
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Fig. 3. — The Tpl distribution of identified sources. A normal 
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each source, the logarithm of V var is expected to follow 
a half- normal distribution. In contrast, the distribution 
of Vvar for compact objects shows a large spread and de- 
pends on the subclasses clearly, as discussed below. 

The distribution of V var for each subclass is given in 
Figure O We see that Orion variables are more variable 
(with m — 2.71) than other stars (such as the main se- 
quence stars, which have m = 1.45). The variable stars, 
which are classified mostly due to their optical variability, 
have m = 1.51, close to that of main sequence stars. The 
three subclasses of Seyfert galaxies (Syln, Syl, and Sy2) 
overall show similar variability, with m = 1.56, 1.57, and 
1.49, respectively. However, some narrow- line Seyfert 1 
galaxies show very large variability (a f actor of 154 for 
PHL 1092; see also iMiniutti et all 120091 LINERs have 
m = 1.16, which is much lower than those of other AGN, 
but we cannot exclude the probability that it is due to 
contamination of stea dy circumnuclear diffuse emi ssion 
in these systems (e.g.. iGonzalez-Martm et aT1l2009l) . 

The large spread of the distribution of V var for com- 
pact objects mostly comes from accreting compact ob- 
jects (BHBs, bursters, accretion-powered X-ray pulsars, 
and CVs; Figures |4][5]) . They were observed to vary by 
factors up to a few thousands. Some of them are known 
to vary much more than measured here (e.g., Aql X-l). 
This can be explained by a selection effect because XMM- 
Newton cannot observe bright sources in imaging modes. 
For isolated NSs, the rotation-powered pulsars and ther- 



mally cooling isolated NSs have m — 1.11 and 1.15, re- 
spectively (their median values of S val are only 0.79 and 
1.35, respectively). PSR J1302-6350 (i.e., PSR B1259-63, 
#116117) is the only one with significantly large vari- 
ability (a factor of 31.2). However, it is a binary system, 
and its X-ray emission mechanism is probably different 
from most rotation-powered pulsars (|Chernvakova et al.l 
2006). Different from the above two classes of isolated 
NSs, the majority of magnetars are observed to be vari- 
able, with some varying by factors of a few tens. 

3.1.4. Multi-wavelength cross- correlation 

We found optical counterparts from the USNO-B1.0 
Catalog for about 75%, 90%, and 56% of the identified 
stars, AGN, compact objects, respectively. The corre- 
sponding percentages of IR counterparts found from the 
2MASS PSC are 100%, 40%, and 56%, respectively. In 
Figure El we show the X-ray-optical and X-ray-IR po- 
sitional separations for the identified stars and AGN. 
About 90% of the optical/IR counterparts have sepa- 
rations <l'/5 from the X-ray positions. Separations are 
expected to follow the Rayleigh distribution if both R.A. 
and Decl. have a constant separation error <7 sep (the 
quadratic sum of the X-ray and optical/IR positional er- 
rors). The value of r7 sep for each Rayleigh distribution 
overplotted on the distribution of the separations in Fig- 
ure [6] (dotted line) was inferred from the median of the 
separations, which is 1.18 times a scp for the Rayleigh dis- 
tribution. We obtained cr se p=0 / . / 42 and 0'.'52 for the opti- 
cal counterparts of AGN and stars, respectively. The cor- 
responding values for the IR counterparts are <7 sep =0'.'42 
and 0'.'51, respectively. The optical/IR counterparts of 
stars on the whole have larger separations from X-ray po- 
sitions than those of AGN, probably due to stellar proper 
motions. The above results indicate a high astrometric 
accuracy of our source sample. 

The identified compact objects for which we found the 
optical/IR counterparts are mostly CVs and accretion- 
powered X-ray pulsars. We found no optical/IR coun- 
terparts for the seven thermally cooling isolated NSs. In 
fact, their optical counterparts are known to be very 
faint , with S-band magnitudes about 26 (jPires et al.l 
2009). Some optical/IR counterparts are found for 
magnetars, but they all have large separations (>1'.'9) 
from the X-ray positions and are mostly in the crowded 
fields of optical/IR sources, thus probably spurious. 
Indeed, magnetars generally have i?-band magnitudes 
larger than 24 and -K~ s -b and magnitudes around 19 
(|Woods fc Thompson! 120061 ) . below the detection limits 
of the USNO-B1.0 Catalog and the 2MASS PSC. We 
found few optical/IR counterparts for the 21 rotation- 
powered pulsars in our sample either. PSR J1302-6350 
(i.e., PSR B1259-63, #116117) and PSR J1023+0038 
(#83189) are the only two confident cases with optical 
and I R counterparts, but both of them are bina ry sys- 
tems (John ston et~atlll992t [Archibald et al.ll2009H . 

Figure [7] plots the X-ray flux versus the optical flux 
for identified sources with optical counterparts. It shows 
that the identified stars span higher optical fluxes than 
AGN and compact objects, but they all have the low- 
est optical fluxes around 10~ 13 5 erg s _1 cm~ 2 , about 
the detection limit in the USNO-B1.0 Catalog. Figure[5] 
plots the X-ray flux versus the IR flux for identified 
sources with IR counterparts. We see that the identi- 
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Fig. 4. — The X-ray flux variation factor versus the maximum flux in 0.2-4.5 keV and the distribution of the flux variation factors. Two 
bursters, #151436 (V var =1.8 x 10 4 ) and #163131 (\/ var =3.0 x 10 4 ), are outside the plotting range in (c). Half- normal distributions of the 
flux variation factor logarithms using the same median as the data are overplotted (dotted lines). 



fied stars also span higher IR fluxes than AGN and com- 
pact objects. Moreover, they all have IR fluxes >10~ 125 
erg s _1 cm -2 , brighter than the 3cr detection limit, while 
many identified AGN and compact objects have IR fluxes 
lower than this limit . 

In the IR color-color diagram in Figure |9] for sources 
with IR counterparts, we see that the identified AGN 
concentrate within the big dashed squ are, which is in 
agree ment with previous studies (e.g., iHvland fc Allen! 
1982). The IR counterparts of the identified stars and 
compact objects span a larger range of colors. Compared 
with the locus of main-sequence stars (the green solid 
line in Figure IH]) and the reddening vector (the arrow), 
their large spread in IR colors is probably due to strong 
extinction in some of them. 

Figure [10] plots the X-ray-to-IR versus X-ray-to-optical 
flux ratios for all identified sources. For sources without 
optical/IR counterparts, their optical/IR flux upper lim- 
its are used and they are circled in this figure. Sources 
with neither optical nor IR counterparts fall onto the 
dotted reference line. The corresponding distributions 
of the X-ray-to-optical and X-ray-to-IR flux ratios are 
shown in Figure [IT] For stars and AGN, we separate 
the sources with (solid line) and without (dotted line) 
optical/IR counterparts (there is no distribution plot for 
stars without IR counterparts because they all have IR 
counterparts). As there are not many compact objects, 
they are not separated but combined to create a sin- 
gle distribution plot. In Figures [TO] [UJ we see that 
stars are generally separated from AGN and compact 
objects in terms of the X-ray-to-IR flux ratio (refer to 
the dashed reference lines; see also Figure [8]). Stars have 
log(^ x /J ; i R )<-0.9, AGN -0.9<log(F x /F m )<2.5, and 
compact objects \og(Fx / Fir)>0 .5 . AGN tend to have 
higher X-ray-to-optical flux ratios (\og(Fx/ Fo)> — 1.0) 
than stars (\og(Fx/ Fo)< — 1.0) for sources with optical 
counterparts. However, many stars, especially Orion 
variables, have no optical counterparts. These stars are 
probably subject to strong extinction so that their op- 
tical counterparts fall below the detection limit of the 
USNO-B1.0 Catalog, making their (apparent) X-ray-to- 
optical flux ratios as large as those of AGN and compact 
objects. 

Figure HT1 shows that, compared with AGN with opti- 



cal counterparts, AGN without optical counterparts on 
the whole have slightly larger X-ray-to-optical flux ratios, 
which were obtained using an (assumed) optical upper 
limit of R2 — 21.0. This is probably due to the selection 
bias that AGN without optical counterparts are far away 
and tend to be included by us if they have high X-ray-to- 
optical flux ratios and are thus bright in X rays, and/or 
the selection bias that AGN with optical counterparts 
include sources with strong star-forming activity in the 
host galaxies and thus with large star light contamina- 
tion. 

3.2. The Source Type Classification Method 

We were able to identify the source types from the 
literature for only about one third of our sources. For 
the rest, we classified them as one of three candidate 
classes: stars, AGN, and compact objects. We derived 
the following source type classification method. 

Sources with log(Fx/i 7 iR)<— 0.9 are classified as stars, 
except those with HR1<0.3 and those in the direction 
of centers of galaxies (from the RC3, the 2MASS XSC, 
the SDSS, etc.). Sources are treated as in the direction 
of galactic centers if the separation is <4". Sources with 
\og(Fx/ Fiji)>— 0.9 are classified as stars if they have X- 
ray flares detected. 

The remaining sources are assumed to be either 
AGN or compact objects. Sources with HR1<— 0.4, 
HR2<-0.5, HR3<-0.7, or HR4<-0.8 (soft sources), 
sources with -0.1<HR3<0.5 and -0.25<HR4<0.1 
(hard sources like accretion-powered X-ray pulsars), 
sources with Vv ar (0.2-4.5 keV)>10, and sources with 
log(i 7 x/-FiR)>2.5 are classified as compact objects and 
denoted as strong candidates. Sources with only a small 
fraction of detections marginally passing the above X- 
ray color selection criteria are not classified as compact 
objects using the colors. To decrease the probability 
of missing compact objects, remaining sources that are 
probably non-nuclear extragalactic sources (a/i?25 < 2.0 
and a > 4") or probably in the Galactic plane (\b\ < 10°) 
are also classified as compact objects but denoted as weak 
candidates. 

Applying the above classification scheme to the iden- 
tified sources we found that it can identify about 
99%, 93%, and 70% (88%) of stars, AGN, and strong- 
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Fig. 5. — The distribution of the X-ray flux variation factors for 
each source type, with the median value also given in parentheses. 
The last bin in each panel includes all sources with variation factors 
>9. 

candidate (plus weak-candidate) compact objects, re- 
spectively. The corresponding probabilities of false iden- 
tification, defined as the ratio of the number of false 
identifications to the total number of identifications, are 
about 1%, 3%, and 18% (32%), respectively. 

3.3. Properties of Candidate Sources 

The results of applying the source type classification 
method outlined in Section[3T2"lto the unidentified sources 
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Fig. 6. — The X-ray-optical and X-ray-IR separations versus the 
maximum 0.2—12.0 keV flux (left panels) and the distributions of 
the separations (right panels) , which are approximated as Rayleigh 
distributions (dotted lines; see text). 



are given in Tables [2] and |4j For very few sources, our 
classification also used specific source properties (such 
as dips) reported in the literature. The relevant refer- 
ences have been included in Table |U We obtained 644 
candidate stars. For candidate AGN, we differentiated 
the 84 coincident with (extended) galaxies (G) from the 
rest (1292). For compact objects, 202 are strong candi- 
dates ("CO" in Tables Hand ID), 320 are candidate non- 
nuclear extragalactic sources (XGSs), and 151 are can- 
didate Galactic-plane sources (GPSs). These sources do 
not include the 16 SNR, the 19 "Mixed" sources and the 
100 candidate ULXs in Tables H and g] Among the 100 
candidate ULXs, 11 were identified in this work. They 
were selected under the conditions of the maximum 0.2- 
12.0 keV luminosity above 10 39 erg s _1 (6 above 2xl0 39 
erg s _1 ), no optical counterparts within 2" of the X-ray 
position from the USNO-B1.0 Catalog or the SDSS, and 
a/R 25 < 1.0 (source #106200 has a/R 25 = 1-99 and 
was included considering its large variability (V var (0.2- 
4.5 keV)=10.9). In Table [H we list the characteristics 
(being soft, highly variable in X-rays, etc.) of the candi- 
date compact objects and ULXs that differentiate them 
from AGN. 

Figures [T2TtT71 show the X-ray color-color diagrams, 
long-term X-ray variability, the X-ray versus optical 
fluxes, the X-ray versus IR fluxes, IR color-color dia- 
grams, and the X-ray-to-IR versus X-ray-to-optical flux 
ratios, respectively, for these candidate sources. They 
correspond to Figuresl2l [4l and [TT-FTOl for identified sources, 
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Fig. 7. — The X-ray (0.2-12.0 keV; maximum) versus optical (USNO-B1.0 i?2-band) fluxes of identified sources with optical counterparts. 
The dotted reference lines, the same for all panels, are plotted to indicate the separation of stars from other source types. 
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Fig. 8. — The X-ray (0.2-12.0 keV; maximum) versus IR (2MASS iC s -band) fluxes of identified sources with IR counterparts. The dotted 
reference lines, the same for all panels, are plotted to indicate the separation of stars from other source types. 
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Fig. 9. — The 2MASS color-color diagram of identified sources. Only the counterparts detected in all three 2MASS bands are plotted. 
One Orion variable (#169259, with J — H = 4.84 and H — K B = 3.51) is outside the plotting range. The big dashed square, which 
encloses the data po ints of mos t AGN , is repeated in all panels for reference. The green solid line represents the locus of main-sequence 
stars, obtained from Koornnccf (1983) a fter c olor transformation (http://www.ipac.caltech.edu/2mass/rclcascs/allsky/doc/sec6_4b.html). 
Reddening by Ay=5 (Rieke &; Lcbofsky 1985) is indicated by an arrow. 
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respectively. The candidate and the identified sources 
generally show similar properties in these plots. There 
are some subtle differences. Comparing Figures H4HT51 
with Figures [7H3 we see that the candidate sources on 
the whole are fainter, mostly with the maxiaml 0.2-12.0 
keV flux <10~ 12 erg s -1 cm~ 2 . This can be explained, 
as brighter X-ray sources tend to have brighter optical 
counterparts, which makes the identification, typically 
through optical measurements, easier. 

Comparing Figures [TU] and [T71 we see that there 
are more candidate stars than identified ones with 
log(Fx/i 7 iR,)>— 0.9. This may be because some of them 
are faint in IR and optical, hense hard to identify, while 
X-ray flares are more effective to identify these sources. 
Sources #173834 and #6007 are two extreme cases. We 
measured log (Fx /-Fir) = 0.63 for source #173834, due to 
a large flare lasting ^20 ks in observation 0200450401. 
It was not detected in five other observations, and we 
obtained a long-term flux variation factor V var (0.2-4.5 
keV)>53. Our classification as a star is also supported 
by appearing within the Cygnus OB2 star forming re- 
gion. Source #6007 has \og{F^/F m ) = 0.65, due to 
a large flare lasting ^8 ks in observation 0505720501. 
It appears within M 31, and it was a lso classifed as a 
foreground star by iStiele et"al"1 (|2008l ) due to a (much 
weaker) flare in an earlier observation. We measured 
K, ar (0.2-4.5 keV)>198. Sources with extreme properties 
warrant future investigation. 

The 84 candidate AGN that are coincident with the 
centers of the extended galaxies on the whole have 
brighter optical and IR counterparts than other candi- 
date AGN (Figures [T41H5)) . They also have lower X- 
ray-to-IR and X-ray-to-optical flux ratios, which separate 
them from other candidate AGN in Figure [TT1 probably 
indicating large contamination of IR and optical counter- 
part fluxes from the starlight. In the X-ray color-color 
diagrams in Figure I12[ some of them are also separated 
from other candidate AGN and are more in the region 
occupied by stars, indicating that their X-ray emission 
may be dominated by the hot gas inside the hos t galaxies, 
as oft en seen in nonactive early-type galaxies (|Fabbianol 
[1981 . 

The 16 SNRs are remarkably similiar to each other. 
They all have X-ray colors similar to stars (Fig- 
ure [T2"j) . Their X-ray emission is probably due to ex- 
panding SNR sh ock heating the sur rounding interstel- 
lar medium (e.g.. Iltoh fc Ma sai 1989!) and can be fitted 
with a bremsstrahhm g model with temperatures <1 keV 
(|Pannuti et al.ll2~000D . They either (for most cases) have 
no IR or optical counterparts found from the USNO-B1.0 
Catalog and the 2MASS PSC, or have counterparts at 
large separations from the X-ray positions (>2"), thus 
probably spurious. They are very steady, with the 0.2- 
4.5 keV flux variation factors all <1.8 (the median is 1.19, 
Figure [5]) and the variability significance S var all <4. 

For the 100 candidate ULXs, most of them are in the 
region occupied by AGN in the X-ray color-color dia- 
grams in Figure[T51 while nearly 30% of them have detec- 
tions entering the soft regions (HR1<— 0.4, HR2<— 0.5, 
HR3<-0.7, or HR4<~0.8) where AGN are hardly seen. 
The U LX sample are often contaminated by background 
AGN (jLopez-Corredoira fc Gutierrez! [2001 . and we ex- 
pect that some of these 100 candidate ULXs are in fact 
AGN. Figure [18] plots the X-ray flux variation factor 



(0.2-4.5 keV) versus the maximum luminosity (0.2-12.0 
keV) and the distribution of the flux variation factor for 
these candidate ULXs. There are 13 candidates (found 
from the literature) with the maximum 0.2-12.0 keV lu- 
minosity slightly less than 10 39 erg s _1 , the lowest limit 
generally used to define ULXs. This discrepancy can 
be explained by various factors, such as source variabil- 
ity and different source distances used in the literature 
from u s (we mos t ly ref er the distances to lLiu fc Bregmanl 
(|2005l) and Eh] (|2011ft ). About 15% of these 100 can- 
didate ULXs have Kar(0.2-4.5 keV)>10, but most of 
them have the maximum 0.2-12.0 keV luminosity lower 
than 4xl0 39 erg s _1 . We obtained a median of 1.73 for 
the 0.2-4.5 keV flux variation factors, which is slightly 
larger than that of AGN but smaller than those of vari- 
ous types of accreting compact objects (Section 13 . 1 .3|) . 
Although we do not have enough BHBs to obtain a 
meaningful value of their average variation factor, we 
expect it to be much larger than that of these candi- 
date ULXs, considering that m ost BHBs are transients 
(jMcClintock fc Remillardll2006l ). 

For the 202 strong candidate compact objects, they 
are spread across the X-ray color-color diagrams (Fig- 
ure [T2]) ; they include some supersoft X-ray sources and 
some hard sources which are probably accretion-powered 
X-ray pulsars. The weak candidate compact objects in 
the Galactic plane are mostly highly absorbed, occupying 
the upper right corner in the HR1-HR2 diagram. 

Several special objects have very low X-ray-to-IR flux 
ratios in the right panel in Figure [T7l though they prob- 
ably contain compact objects. The four candidate sym- 
biotic stars, i.e., AG Draconis (#143273), 4U 1700+24 
(#152155), Z And (#189311), and SMC Symbiotic Star 
3 (#194075), all have log(F x /i ; iR)<-1.0 (they are in- 
cluded in the "Mixed" class). One supergiant fast X- 
ray transient (SFXT), i.e., IGR J17544-2619 (#158580), 
has log(Jx/-FiR,)=— 1.0 (it is included in the "CO" class, 
iGonzalez-Riestra et al.ll2004D . There is another SFXT 
IGR J11215-5952 (#206717), which has been included in 
the class of accretion-powere d X-ray pulsars d ue to the 
detection of the spin period (jSidoli et al.l 120071 ) and has 
log(Fx/-FiR,) only —0.1. We also detected flares lasting 
~1 h from these two SFXTs. Their flares, however, show 
complicated structure, in contrast with the typical fast 
rise and slow decay of stellar X-ray flares. Observations 
of very hard spectra, short pulsation periods and extreme 
long- term variability as in IGR J11215-5952 (Sido li et al.l 
2007) can help to differentiate them from stars. 

4. DISCUSSION AND CONCLUSIONS 

We have systematically studied the properties of 4330 
point sources from the 2XMMi-DR3 catalog. These 
sources represent a high-quality X-ray source sample, 
each with multiple observations and at least one detec- 
tion with S/N>20. They have very high astrometric ac- 
curacy (99% with the positional error <1"). For about 
one third of them we have obtained reliable source types 
from the literature. They correspond to various types of 
stars, AGN and compact object systems containing WDs, 
NSs, and stellar-mass BHs. We have studied their prop- 
erties in terms of the X-ray spectral shape, X-ray vari- 
ability, and multi-wavelength cross-correlation. We find 
that 99% of stars can be separated from other sources us- 
ing the X-ray-to-IR flux ratio and the X-ray flare. Stars 
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Fig. 10. — The X-ray-to-IR versus X-ray-to-optical flux ratios of identified sources. The dotted and dashed lines are plotted for reference. 
Data points circled have no optical (to the right of the dotted line) and/or IR (to the left of the dotted line) counterparts found, and their 
optical/IR fluxes were calculated using detection limits (see text). 

also occupy well defined regions in the X-ray and IR 
color-color diagrams where sources of other types, espe- 
cially AGN, seldom visit. Comparing AGN and compact 
objects, the X-ray spectra of AGN are remarkably similar 
to each other, with Tpl concentrated around 1.91±0.31, 
and their long-term flux variation factors have a median 
of 1.48 and 98.5% less than 10, while 70% of compact 
object systems can be very soft or hard, highly variable 
in X-rays, and/or have very large X-ray-to-IR flux ratios, 
separating them from AGN. Using the above results, we 
have derived a source type classification method to clas- 
sify the rest of the sources in our sample. 

In our source type classification method, we select 
sources with log(Fx/-FiR)<— 0.9 and/or X-ray flares to 
be stars unless they have HR1<0.3 and/or are coinci- 
dent with the centers of extended optical/IR sources, 
presumably galaxies. Due to the requirement that the 
source has at least one detection with S/N>20 in our 
source selection criteria, our sources have maximum 0.2- 
12.0 keV flux 10 -14 to 10 -9 erg s _1 cm~ 2 , corresponding 
to relativelly bright sources in the 2XMMi-DR3 catalog. 
For too faint sources when their IR counterparts fall well 
below the 2MASS detection limit, it becomes hard to es- 
timate their log(Fx / -Fir) ■ For them, the X-ray color and 
the flare (about 20% of stars have flares detected by us) 
will be more effective to separate stars from other types 
of sources than the X-ray-to-IR flux ratio. To separate 
compact object systems from AGN, we rely on X-ray col- 
ors, X-ray long-term variability and the X-ray-to-IR flux 
ratio. There are many other strong indicators of com- 
pact object systems, such as dips, eclipses, pulsations 
(coherent or noncoherent) and bursts, but they normally 
require high-quality data and detailed analysis. 

We find 644 candidate stars, 1376 candidate AGN and 
202 candidate compact object systems, whose spurious 
probabilities arc estimated to be about 1%, 3% and 
18%, respectively, based on the test of our classifica- 
tion method on the identified sources. There are still 
320 appearing to be associated with nearby galaxies and 
151 appearing in the Galactic plane, which are probably 
compact object systems or background AGN. Our source 
sample also contain 100 candidate ULXs, which include 
11 new ones. With long-term flux variation factors hav- 
ing a median of only 1.73 and 85% less than 10, they are 
much less variable than other accreting compact objects. 




-2 2 

-ray to optical/infrared flux ratio 



Fig. 11. — The distributions of the X-ray-to-optical (top three 
panels) and X-ray-to-IR flux ratios (bottom three panels) of iden- 
tified sources. For stars and AGN, we separate the sources with 
(solid line) and without (dotted line) optical/IR counterparts, but 
not for compact objects. The dashed lines in the distribution plots 
of X-ray-to-IR flux ratios are for reference. 
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Fig. 12. — The X-ray color-color diagrams, the same as Figure[2]but for candidate sources. 
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Fig. 13. — The X-ray flux variation factor versus the maximum flux in 0.2—4.5 keV and the distribution of the flux variation factor, the 
same as Figure [4] but for candidate sources. 
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Fig. 14. — The X-ray (0.2-12.0 keV; maximum) versus optical (USNO-B1.0 i?,2-band) fluxes, the same as Figure [7] but for candidate 
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Fig. 15. — The X-ray (0.2-12.0 keV; maximum) versus IR (2MASS i\s-band) fluxes, the same as Figure|8]but for candidate sources. 
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Fig. 16. — The 2MASS color-color diagram, the same as Figure [9] but for candidate sources. Two stars, #146483 (J — H = 4.92, 
H — K s = 3.14) and #146372 (J - H = 5.74, H — K s = 3.91), are outside the plotting range. 
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Fig. 17. — The X-ray-to-IR versus X-ray-to-optical flux ratios, the same as Figure ITOl but for candidate sources. 
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Fig. 18. — The X-ray flux variation factor (0.2-4.5 keV) versus 
the maximum luminosity (0.2-12.0 keV) and the distribution of 
the flux variation factor for candidate ULXs. 

In this work, we have also found a few tens of sources 
that show interesting properties but are poorly studied 
in the literature. We will present our detailed study of 
these sources in a companion paper. 

There are other studies using multi-wavelength 
cross-correlation to cla s sify X -ray sources (See, e.g., 

lHaakonsen fe Rutledgel 120091 a nd the references 

therein) . lHaakonsen fc Rutledgel (|2009l ) studied the 



cross-correlation between the RASS Bright Source 
Catalog (RASS BSC) and the 2MASS PSC and showed 
that stars are well separated from other types of sources 
in terms of the X-ray-to-IR flux ratio and the IR color 
J — K s . As sources in the RASS BSC mostl y have the 
0.1-2 .4 keV flux >10~ 12 erg s" 1 cm" 2 (jVoges et all 
1999), our sample extends the fluxes by two orders of 
magnitude below theirs. Although we also find that 
stars generally have lower X-ray-to-IR flux ratios than 
other types of sources, their result that stars have 
J — K s < 1.1 is not well supported by our source 
sample. Our identified and candidate stars have J — K s 
spanning from about 0.0 to 8.0, and about 30% and 
20% of them have J — K s > 1.1, respectively. The high 
value of the J — K s color for these stars is probably 
caused by large extinction. The lack of these highly 
absorbed sources in the RASS BSC is easy to explain, 
as the energy pass of ROSAT is in the very soft X-rays 
(0.1-2.4 keV). Because of these highly absorbed sources, 
we choose to use the K s -han& flux to calculate the 
X -ray-to-IR flux ratio ins t ead o f the J s -band flux used 
in lHaakonsen fc Rutledgel ((2009) . This is also why stars 
in our sample are not well separated from other types 
of sources in terms of the X-ray-to-optical flux ratio, 
thoug h this is observed to be the case in many studies 
(e.g.. lVoges et al.lll999l ). 
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APPENDIX 

ESTIMATE OF THE SYSTEMATIC ERRORS OF X-RAY FLUXES AND HARDNESS RATIOS 

The 2XMMLDR3 catalog only gives the statistical errors for the X-ray fluxes. Their systematic errors are estimated 
as follows. We first considered the systematic difference of the fluxes between pn and MOS. The MOS flux is the 
average of the fluxes of MOS1 and MOS2 (they have very similar responses) weighted by the errors. Using detections 
with both the pn and MOS fluxes and at least one of them above 5 a from our source sample, we determined the 
systematic errors that make 31.7% (corresponding to l-a difference) of detections with the pn and MOS flux difference 
larger than the total error, i.e., 

\F pn - FmosI > [(fpn x F pn ) 2 + (r MO s x F MO s) 2 + (ct p „) 2 + (o-mos) 2 ] 1/2 , (Al) 

where F pn and -Fmos are the pn and MOS fluxes, with the l-a statistical errors a pn and ctmos and fractional systematic 
errors r pn and r pn , respectively. Assuming r pn = tmos, we estimated r pn = tmos == 0.052 and 0.076 for energy bands 
14 and 8, respectively. There are also systematic errors due to different pointing offsets, data modes, etc. for different 
observations and the approximation of the PL spectrum, and they are hard to estimate. In the end, we just assumed 
the fractional systematic errors of the EPIC fluxes in energy bands 14 and 8 to be r = \/2r pn = 0.074 and 0.107, 
respectively. 

The systematic differences of the hardness ratios between the pn and MOS can be estimated in a similar way. Using 
detections with both the pn and MOS hardness ratios and at least one of their statistical errors less than 0.1, wc 
determined the systematic errors that make 31.7% of detections with the difference of the pn and MOS hardness ratios 
larger than the total error, i.e., 

|HRi pn - HRiMOsI > [(Si.pn) 2 + (s*,MOs) 2 + Kpnf + (ff^MOS )*] , (A2) 

where HRi pn and HRzmos are the hardness ratios defined as in Equation Q] using the pn and MOS fluxes, respectively, 
with the statistical errors <7j iPn and ct^mos and the absolute systematic errors Si tPn and mos, respectively. Assuming 
Si, pn = Si, MOS: we obtained Si. pn — s^mos = 0.036, 0.041, 0.014, and 0.030 for i=l to 4, respectively. In comparison, 
the hardness ratios in our definition assume values between —1 and 1. 

SEARCH FOR STELLAR X-RAY FLARES 

We searched for stellar X-ray flares as follows. To minimize the risk of false detection of flares due to flaring 
background, we excluded data in the period of flaring background in the search. We first located the data point i 
at time t{ which has the maximum weighted count rate Ri calculated from its nearby three data points i, i + 1, and 
i + 2. The weighted error of Ri is denoted as Ri^ri- We then calculated the weighted background count rate bi and 
weighted standard deviation bi iCrl using data in the background time windows [U — 40ks, U — t w - m ] and [U + 2t w - m , 
ti + 40ks], where i w i n is equal to 10 ks or one fifth of the observation duration for observations lasting <50 ks. We 
identified a flare candidate at U if Ri — 2Ri C „ > bi + 2&i jC1T and Ri/bi > 5. This process continued, but using only 
data 40 ks outside the previously identified flare candidates. The test on identified stars (Section 12. 4|) showed that 
the above procedure can successfully identify most flares. Some flares (about 10%) are missed mostly because they 
are very weak or the observations are short. A few non-stellar flares from sources with large fast variability (such as 
periodic oscillations in AM Hcrculis objects) were also picked up. Their profiles are generally very complicated and 
very different from the typical profiles of smooth fast rise and slow decay of stellar flares, and we discarded them based 
on visual inspection. In Table |4] we give the maximum ratio Ri/bi and the corresponding observation for sources with 
stellar flares found. For sources that have no flares found using the above search procedure but have some flares from 
our visual inspection, this ratio is set to be one. 
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TABLE 3 

Source Detections and Simple Fits with an Absorbed PL 



SRCID 
(1) 


OBSID SRCIDDet 
(2) (3) 


Rl Rl.crr R2 R2_crr R3 R3_crr R4 R4_crr R5 R5_crr 

(4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 


X 2 dof Nh 
(14) (15) (16) 


»H,le 
(17) 


™H,u< 
(18) 


r PL 

(19) 


r PL,le 
(20) 


r PL ,n< 

(21) 


F PL F PL,le F PL,uc 
(22) (23) (24) 


837 


0101040101 


837 


1.67o-02 7.46c-04 1.40e-02 6.82e-04 7.56e-03 5. 19e-04 2.49e-03 3.11e-04 4.99e-04 2.89e-04 


1.4 


6 0.00 


0.00 


0.03 


2.92 


2.79 


3.23 


3.90e-13 3.03e-134.37e-13 


837 


0306870101 


837 


1.03e-02 3.10e-04 9.23e-03 3.08e-04 5.68e-03 2.58e-041.92e-03 1.60e-04 5.22e-04 1.63e-04 


2.0 


6 0.00 


0.00 


0.02 


2.72 


2.61 


2.95 


2.54o-13 2.07o-13 2.85e-13 


837 


0510010701 


837 


8.55e-03 4.93e-04 7.10e-03 4.96e-04 4.42e-03 4.83e-041.58e-03 3.18e-041.01e-04 2.31e-04 


0.9 


6 0.00 


0.00 


0.03 


2.69 


2.50 


3.13 


1.69c-13 1.35o-13 1.94e-13 


938 


0101040101 


938 


1.78o-03 1.37e-042.23e-03 1.57e-042.55e-03 1.90e-049.24e-04 1.27e-04 5.51e-04 1.33e-04 


1.6 


10 0.00 


0.00 


0.03 


2.04 


1.90 


2.33 


6.75c-14 5.74o-14 7.93e-14 


938 


0306870101 


938 


1.12c-03 1.05e-04 1.92e-031.39e-04 2.26e-031.57e-041.44e-03 1.30e-04 4.61e-04 1.32e-04 


0.9 


6 0.00 


0.00 


0.04 


1.72 


1.58 


1.98 


8.64e-14 7.02c-14 1.01e-13 


938 


0510010701 


938 


1.03c-03 1.72c-04 1.14c-03 1.99c-04 9.62c-04 2.44c-04 3.32c-04 1. 72c- 04 1.51c-04 1.81c-04 


0.6 


6 0.01 


0.00 


0.12 


2.42 


1.96 


3.58 


2.70e-14 1.70e-143.81e-14 


Note 


— This table 


is publi 


shed in its entirety in the electronic edition of the Astrohysical Journal. A portion is 


sho 






regar 




form £ 


nd content. Columns: (1) 



the unique source number that we used to refer the source, (2): observation ID, (3): the unique source number in the 2XMMi-DR3 catalog: (4)-(13): the MOSl-medium-filter equivalent 
counts rates in five basic energy bands and the corresponding errors (in units of cts/s), (14): the reduced % 2 of the fit with an absorbed PL, (15) degrees of freedom, (16)— (18): the 
best-fitting Njj and lower and upper bounds (at a 90% confidence level, in units of 10 22 cm -2 ), (19)-(21): the best-fitting PL photon index and lower and upper bounds, (22)— (24): 
the 0.2-12.0 keV flux from the PL fit and the lower and upper bounds {in units of erg cm -2 s — For a unique source, the detections arc specified by OBSID and SRCIDDet, which 
are "OBS.ID" and "SRCID" in the 2XMMi-DR3 catalog, respectively. SRCID is different from SRCIDDet only when several "unique" sources in the 2XMMi-DR3 catalog arc rematched 
together (see Tabic 1). When the source has no detection in the 2XMMi-DR3 catalog for an observation, the corrcsonding SRCIDDet is empty, and the fit was not carried out. There 
are 17 detections not used in our study (thus not in this table too), because of large discrepancy between pn and MOS measurements or large extended radius used in the 2XMMi- DR3 
catalog (causing fake large fluxes). 



TABLE 4 

General Source Properties and Source Classification 



SRCID 2XMMUDR3 RAdcg 


DEdeg 


NObs 


DObs 


NDet 


DDet 


Type 


TpCat 


Ref 


Ref Type Sr 


cChar a 


Rxo 


Rxir 


Vvarl4 


Dxo 


Dxir 


RC3Name 


RC3SepR 


(1) (2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


8085 2XMM J005412. 9-373309 13.55411 


-37.55272 


4 


1794.51 


4 


1794.51 


CO 


B 






SVE 


1.26 


0.56 


36.33 






NGC300 


1.11 


100081 2XMM J120429. 7+201858 181.12376 


20.31638 


3 


178.09 


3 


178.09 


Sy2 


A 


VV10 


A/S2 




-2.06 


0.11 


2.04 


0.6 


0.8 






101335 2XMM J121028. 9+391748 182.62064 


39.29689 


8 


2169.02 


8 


2169.02 


CO 


B 






S 


0.70 


0.00 


2.63 










107102 2XMM J123103. 2 + 110648 187.76337 


11.11348 


3 


887.86 


3 


887.86 


CO 


B 






S 


0.31 


0.12 


2.11 


3.6 








112394 2XMM J125048. 6+410743 192.70256 


41.12836 


3 


1649.30 


2 


1615.12 


ULX 


B 






VE 


1.67 


0.97 


25.85 






NGC4736 


0.17 


250424 2XMMi .1170058.4-461107 255.24352 


-46.18554 


3 


209.32 


3 


209.32 


Bstr 


A 


2009Ap.I...699. 


.60L Bstr 




1.92 


1.22 


3.15 











NOTE. This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. Columns 

arc as follows. (1): 2XMM-DR3 Unique source index, (2): 2XMMi-DR3 Source designation, (3): Mean source ICRS right ascension (J2000), (4): Mean source ICRS declination (J2000), 
(5)-(6): Number and time span (in units of days, using the observation start time) of XMM-Newton observations, (7)-(8): Number and time span (in units of days, using the observation 
start time) of detections from the 2XMMi-DR3 catalog, (9): Source type, (10) The source type category ( A=identified, B — candidate) , (11)-(12): Reference and the spectral type in 
the reference, (13): The source characteristics based on which we classify the source as a compact object system (instead of star or AGN), (14): X-ray-to-optical flux ratio logarithm 
log(F x /FQ), (15): X-ray-to-IR flux ratio logarithm log( F^ / Fj j^) , (16): Flux variation factor in the 0.2-4.5 keV band, (17): USNO-B1.0 counterpart separation from the X-ray position 
(in units of arcsec), (18): 2MASS counterpart separation from the X-ray position (in units of arcsec), (19): The source name of the RC3 match, (20): The ratio of the angular separation 
to the -D25 isophotc elliptical radius. 

a "S": soft (HRK-0.4, HR2<-0.5, HR3<-0.7, or HR4<-0.8); "H" : hard (-0.1<HR3<0.5 and - . 25 < HR4 < 0. 1 ) ; "V": highly variable V V ar(0.2-4.5 keV)>10; "R" : high X-ray-to-IR 
flux ratio (log( F x / Fj R ) > 2 . 5) ; L 'E" : non-nuclear extra- galactic; "G" : in the Galactic plane (|fa|<10); "L" : literature, which includes more source properties such as dips. 



